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1, INTRODUCTION 

Nowadays, the AC machines have replaced the DC machines in industrial applications because 
of their advantages, such as, the reliability and the lack of commutator and brushes which make them able 
to work under unfriendly conditions. The most popular AC machines are the induction motors (IM) 
and the permanent magnet synchronous motors. They are used in various industrial applications, electric 
vehicles, tools and drives etc. The squirrel cage induction motor in particular, is widely used due to its 
reduced cost and lower maintenance requirement [1]. 

In the early decades, the induction motors have been operated directly from the grid under a fixed 
frequency-speed. Later, with the development of modern semiconductor devices and power electronic 
converters, these machines had become able to operate with adjustable frequency-speed by supplying them 
through a power converter like the voltage source inverter. The employment of the variable speed motor 
drive in open loop may offer a satisfied performance at steady state without need of speed regulation 
for simple applications [2]. But, in cases where the drive requires fast dynamic response and accurate speed, 
the open loop control becomes unsatisfactory. Therefore, it is necessary to operate the motor in a closed loop 
mode. Several techniques have been proposed for this purpose. They are classified mainly into scalar 
and vector controls [3]. 

The vector control, which is known also by the field-oriented control (FOC), was developed 
to overcome the limitation of the scalar control. It was presented in the 1970s by Hasse [4] and Blaschke [5] 
to provide an independent control of torque and flux in similar way to the separate excitation DC machine. 
The vector representation of the motor quantities makes it valid to work in both steady and dynamic 
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conditions, this achieves a good transient response. In the control algorithm of FOC _ based 
on the transformation to the synchronous frame, all quantities will appear as DC quantities. 

However, the standard FOC method suffers from high phase current, line voltage and torque ripples. 
Many modified FOC schemes have been proposed in order to reduce ripples. The insertion of the space 
vector modulation (SVM) was a very useful solution [6]. This method, known by SVM-FOC, reduces 
the high ripples level in spite of its complexity. SVM-FOC algorithm uses linear PI torque and flux 
controllers to generate the reference control voltages [7]. 

The initial mathematical model based on conventional analytic methods usually contains 
approximate hypothesis and unmodeled dynamics. Moreover, it may be affected also by the variation in 
parameters due to environment conditions and the external disturbance during operating. Consequently, 
the use of linear methods such the PI controllers cannot achieve high promising performance. The nonlinear 
controllers can offer several advantages compared to linear control schemes. Among the interest researches 
in the field of nonlinear control techniques are the sliding mode control (SMC) [8] and the backstepping 
control (BSC) with integral action [9]. 

Moreover, the sensorless control is another major issue in control domain. The developed control 
schemes using advanced strategies such as the SVM and nonlinear techniques require an accurate speed and 
flux measurement or estimation for closed loop control design [10]. The use of sensors has several downsides 
like high cost, fragility and low reliability. Furthermore, the physical environment sometimes, does not allow 
to use sensors. Due to the multiple variables and nonlinearity of induction motor dynamics, the estimation of 
the rotor speed and flux without the measurement is still very challenging subject. Various sensorless 
approaches have been proposed in the literature [11]. 


2. INDOCTION MOTOR MODELING 

Developing new control laws of observation and optimization for the IM naturally imposes not to 
neglect the modelling aspect of this machine. The modelling of AC machines is based mainly on the work of 
G. Kron, who gave birth to the concept of generalized machine [12]. Park’s model is a special case of this 
concept. It is often used for the synthesis of control laws and estimators. Described by a non linear algebra- 
differential system, Park’s model reflects the dynamic behaviour of the electrical and electromagnetic modes 
of the asynchronous machine. It admits several classes of state representations. These model classes depend 
directly on the control objectives (torque, speed, position), the nature of the power source of the work 
repository and the choice of state vector components (flux or currents, stator or rotor). The IM drive can be 
described by the following state equations in the rotating field reference frame. 


X = AX + BU 


Y = CX ) 


Where X, U and Y are the state vector, the input vector and the output vector respectively. With: 


X= ligsigsParParl ;U= [Wastgs] ;¥ = liasigs] 


1 





oLs 
A —-Ws —a, —a3.Wr a2 . B=] 0 iy] | C= be 0 0 "| 
a, =O —as Wr |’ ols} ’ 01 0 0 
0 ai —W, —ac 0 0 
0 0 
mda, joe, Se, = eS eS = 
Gli.  Ghdly olgLeTy oLsLy tc i Ry 


3. DIRECT FIELD ORIENTED CONTROL BASED BACKSTEPPING APPROACH 
3.1. Main objective of vector control 

The main objective of the field-oriented control is to set a simple model of the IM drive in which 
the flux component and the current magnitude are separated. For a DC machine, the field inductor 
and the induced current are naturally decoupled and orthogonal. In the case of an IM, it is difficult 
to distinguish the torque-producing current and the flux-producing current, which are strongly 
coupled [13, 14].The idea of the oriented field is to choose a two-phase rotating axis system (d, q) mounted 
on the rotor, stator or air gap, allowing decoupling the torque and the flux. Rotor flux orientation control 
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is the most widely used because it eliminates the influence of rotor and stator leakage reactance and gives 
better results than methods based on stator flux or gap orientation [11]. The condition of this alignment 
is expressed by the following relation: 


Dar = Or ; Dor = 


The direct version of the vector control requires a good knowledge of the module and the phase 
of the flux. The first idea is to access the flux from sensors located in the air gap of the machine. 
However, the introduction of these sensors causes an increase in volume and weakening of the machine [15]. 
In addition, these sensors are sensitive to mechanical shocks and temperature variations. Instead, measured 
quantities such as stator voltages and currents can be used to estimate or observe the rotor flux. In this type 
of control, the norm of the flux and its position used for the transformation of the coordinates 
are determined as follows: 





by = Vix + Pfr) (2) 
= er 
6, = Arctg( ~) (3) 
With: 
Ly 
Dor a Ts J (Ves _ Rgigs)dt ~ OLsigs) 
Ly 
Der = i J (Ves _ Reig, )dt = OLgigs) 


Note that the essential problem of the vector control is to determine the position and the norm of 
the rotor flux. 


3.2. Backkstepping approach applied to direct vector control of IM drive 

The backstepping approach is a technique of systematic and recursive synthesis intended 
for the class of nonlinear systems, the basic idea behind backstepping is to recursively select some 
appropriate state functions as virtual command inputs for first-order subsystems of the global system. In each 
step, an extended Lyapunov function is associated to ensure the stability of the entire system [16]. 
Step 1: 

In this step, the speed and the rotor flux are regulation components; we define the regulation errors 
e, and ez by: 


ea = OF- 0 
ep = Oar — Par 


Since the objective requires that both errors converge to zero, we choose igs and igs as virtual 
commands for the system, for this, we consider the candidate Lyapunov function: 


1 

y= 5 [ej + e3] (4) 
Its derivative is: 
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Where k, and kg are positive constants. 


So that the derivative of the Lyapunov function is negative, the virtual commands representing 
the stabilizing functions can be chosen as: 


i; = —— [ke + 0" + 


f 
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wk Tr . 1 
las = Lm [kze, + Dar oF 7, Parl (7) 


We then obtain: 


V; = —k,e? = k,e% <0 


The virtual commands in (6) and (7) are chosen to satisfy the pursuit objectives and are also considered 
as references for the next step. 
Step 2: 

Now the currents ig, and ig, are the new control components considered as virtual commands 
of this step, we define the regulation errors e3 and e,: 


* 


= . . eo ke axk . 
Therefore, error dynamics e3 and e, can be expressed by: 


; : 1 


Where: 
; ; Lm ldsiqs 
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Consider the following Lyapunov candidate function: 
v= = [e3 + e% + e3 + 4] (10) 
Its time derivative is given by: 


° - 1 * 
V, = —k,e? — k,e3 — k3e% — ky,e4 + €3 | 3e; + las — Wy — = Vas| + e,[kKye, +14, — V2 — 
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Where kg and Ky are positive constants. Then, the following control laws are deduced: 
Vgs = OLs[K3e3 + 15 — Wi] (11) 
Vag = OLs[K4e, +15, — We] (12) 
Which makes the derivative of the extended Lyapunov function Vznegative: 
V, = —k,e? — k,e2 — k3e2 — kye% < 0 


Proposition 1:Consider the model of the IM (1). If the speed and flux controllers synthesized 
by the backstepping method are given by (11) and (12) respectively, then these commands are such that 
the flux and speed converge asymptotically to their desired values. 


3.3. Backstepping approach with integral action for robust control 

It is clear that the structure of the controller generated by the classic version of the backstepping 
is composed of a proportional action, to which is added a derivate action on the error. Such a structure makes 
the system sensitive to measurement noises [17]. The absence of the integrator action leads to the appearance 
of a non-zero constant static error. The solution of this problem is the design of a new version 
of the backstepping with an integral action. Integrators should be integrated into the machine model 
and apply the conventional backstepping method to this new model. The integral action will be automatically 
transferred from the model to the control law [18]. 
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Speed tracking error is defined by: 

Z, =e, + §,f e,dt 
With 8, is a positive constant and 6, { e,dt is the integral action added to the backstepping control to ensure 


convergence of the speed tracking error to zero. 
Consider the candidate function of Lyapunov and its derivative given by: 


1 
Vy = 524 (13) 


V; S=kize +z; [sea + 0* — UDarigs + “ + 0 + 81e;| 


According to the backstepping method and in order to ensure the speed tracking stability, the virtual control 
igs 1S given by the following equation: 


ag 1 Sosa 
GS ae [k,z, + 0% + a + 5,e,] (14) 


We then get: 


V; = —k,z? < 0 


The state ig, is then used as an intermediate command to ensure stability. The time derivative of 2; 
is given by: 


Zz = las — lgs 
Consider the extended Lyapunov candidate function and its derivative given by: 


V2 =V, +23 (15) 


: 1 
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By choosing the control law Vg, as: 


Vgs = OLs[K3z3 + 145 — Wi] (16) 
We find that: 
V, = —k,z? —k3z2 < 0 (17) 


4. SPACE VECTOR MODULATION ALGORITHM FOR RIPPLES REDUCTION 

In this FOC strategy, the pulse width modulation (PWM) is replaced by the space vector modulation 
(SVM) for the purpose of voltage vector selection. This control schemes can preserve a constant switching 
frequency, consequently it reduces the high torque and phase current ripples which is the main drawback 
of the conventional FOC. SVM was first presented in the second half of the 1980s as an alternative strategy 
to the basic Pulse Width Modulator (PWM). Since then, a lot of work has been done on its theory and 
implementation. SVM bases on the space vector representation of the voltage inverter output. It has several 
advantages such as minimizing ripple and Total Harmonic Distortion (THD) and switching losses [19]. 

SVM is different from the conventional pulse width modulation (PWM). It relies on the space vector 
representation of the inverter output. There are no separate modulators for each phase. The reference voltages 
are given by space voltage vector (i.e. voltage vector components in the complex plan). The principle 
of SVM is the prediction of inverter voltage vector by the projection of the reference vector Y% between 
adjacent vectors corresponding to two non-zero switching states [20]. For two-levels inverter, the switching 
vectors diagram forms a hexagon divided into six sectors, each one is expanded by 60° as shown in Figure 1. 
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The application time for each vector can be obtained by vector calculations and the rest of the time period 
will be spent by applying the null vector. When the reference voltage is in sector | as shown in Figure 2, 
the reference voltage can be synthesized by using the vectors V,, V2, and Vo(zero vector) [21]. 
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Figure 1. Diagram of voltage space vector Figure 2. Reference vector as a combination 
of adjacent vectors at sector 1 
The volt-second principle for sector 1 can be expressed by: 
V3 T, = V,T, + V2T2 + VoTo (18) 
de — Ty + T, + To (19) 


T,, T, and Tpare the corresponding application times of the voltage vectors respectively. Tis the sampling 
time. The determination of times T,and ‘T,corresponding to voltage vectors are obtained by 
simple projections: 








Tz * * 

T=52 (VOVig — V2Veu) (20) 
V2T 2 yy 

Ip Vee (21) 


With Vj. is the DC bus voltage.The calculation of switching times (duty cycles) is expressed 
as follows [22, 23]. as shown in Figure 3 andsummarizes switching times (output) for each sector 
as shown in Table 1. 
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Figure 3. Switching times of sector 1 
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Table 1. Switching times for each sector 


_ Sector ot 2B 
Sa Thon Taon Taon Toon Thon Toon 
Sp Taon Toon Thon Thon aa Taon 
Se T, Ty T, T, Ta Ty 





5. ADAPTIVE FUZZY LUENBERGER OBSERVER FOR FLUX AND SPEED ESTIMATION 

The Adaptive flux observer is a deterministic type of observers based on a deterministic model 
of the system. In this work, the adaptive Luenberger state observer (ALO) is used to estimate the flux 
components and rotor speed of the IM drive by including an adaptive mechanism based on the Lyapunov 
theory [24] and fuzzy regulation. In general, the equations of the ALO can be expressed as follow: 


ts =AX+BU+L(¥-¥) (25) 
Y=C.X 
The symbol * denotes the estimated value and L is the observer gain matrix. The mechanism 


of adaptation speed is deduced by Lyapunov theory. The estimation error of the stator current and rotor flux, 
which is the difference between the observer and the model of the motor [25], is given by: 


é=(A-—L.C).e+AA.X (26) 
0 0 0 K. Aw, 
_ a _ ~ |0 0 —-K.Aw, 0 
Where e = X — X AA =A-A= 0 0 0 Nes 
0 0 Aw, 0 


We consider the following Lyapunov function: 


(Aw)? 


V=e.e+ (27) 


Where A is a positive coefficient, its derivative is given as: 
V = e'[(A— LC) + (A— LO)Je — 2KAw,. (Cins Per — Cigs Par) 


With @, is the estimated rotor speed. The adaptation law for the estimation of the rotor speed can be deduced 
by the equality between the second and third terms of (9): 


Or = [2X K. (Cias Per = CigsPar)- dt (28) 


The feedback gain matrix L is chosen to ensure the fast and robust dynamic performance of the closed loop 
observer. 

Lh 
Lok 
i aly (29) 


I, |; 


With 1,, I,, lz and 1, are given by: 


2 

h=(i-D0+D 3 b=-Wu-D.G ; b=SX(y-K2)+S2 +5 | 

x, NY oe 

4 ie oor 
k, is a positive coefficient obtained by pole placement approach; a wise choice was made for its value 
in order to guarantee a fast response. Fuzzy logic control (FLC) in adaptation mechanism replaces 
conventional control and it gives robust performance against parameter variation and machine 
saturation [26, 27]. Figure 4 shows the overall structure of the proposed sensorless robust control [28]. 
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Figure 4. The overall structure of the proposed control 


6. RESULTS AND DISCUSSION 

To demonstrate the effectiveness of the proposed strategy, the control method and the state observer 
are implemented in MATLAB/Simulink. The used induction motor rated power and parameters 
are given in Table 2. 


Table 2. Induction motor rated power and parameters 


Variable Symbole Unit 
Rated power P 3 kW 
Stator resistance Rg 2.20 QO 
Rotor resistance R, 2.68 Q 
Stator inductance Lg 0.229 H 
Rotor inductance Ly 0.229 H 
Mutual inductance Lia 0.217 H 
Inertia moment J 0.047 kg. m? 
Friction coefficient F,. 0.004 kg. m/s 
Poles pair p 2 


6.1. Starting up and steady state with load application 

This section presents the starting up state of the induction motor according to speedstep reference 
of 100 rad/s. Then, a load of 20N.m at t=0.5s is applied and removed at t=0.9s. Figures 5-11 show 
respectively rotor speed response, rotor speed error, rotor flux, rotor flux error, electromagnetic torque, stator 
phase current and line voltage.The figures are specified: (a) for PWM-FOC, and (b) for SVM-FOC. 
The Figures 5-11 show that both techniques prove good dynamic at starting up. We can notice that the speed 
regulation loop rejects the applied load disturbance quickly. The SVM-FOC in Figure 5(b)kept the same fast 
speed response of FOC strategy. Since the same backstepping controller is used for both schemes, there is no 
difference in the transient response. 

Then, Figure 9 illustrates the torque responses with load application. It shows that at the beginning 
the backstepping speed controller operates the system at the physical limit. The PWM-FOC in Figure 10(a) 
shows a chopped sinusoid waveform of current which indicates the high harmonics level, while SVM-FOC in 
Figure 10(b) shows a smoother sinusoid waveform. The speed error between the observed and the motor 
speed does not exceed 3 rad/s in SVM-FOC while it exceeds that in PWM-FOC. The speed response 
shows a fast and accurate tracking of the imposed high-speed trajectory with minimal steady state error 
and rise time. Finally, Figure 1lshows the comparison of inverter live voltage for both techniques. 
For the PWM-FOC, it is observed that the voltage amplitude is higher than SVM-FOC, and this can increase 
the commutation losses. 
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Figure 5. Rotor speed response 
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Figure 8. Rotor flux error 
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Figure 10. Stator phase current 
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Figure 11. Line voltage 


6.2. THD improvement 

This section presents the THD analysis for the electromagnetic torque, the stator phase current 
and the line voltage. The figures are specified: (a) for PWM-FOC, and (b) for SVM-FOC. The torque THD 
is reduced from 64.21% to 54.03%, the phase current THD is reduced from 30.94% to 17.87% and the line 
voltage THD is reduced from 75.63% to 65.54%. The switching frequency of the SVM-FOC technique 
is constant, the reason is that each inverter’s interrupter has a rest moment (switching off) which can reduce 
the commutation losses even in steady state or in transient operation. Figure 12 shows electromagnetic torque 
THD. Figure 13 shows phase current THD. Figure 14 shows line voltage THD. 
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Figure 12. Electromagnetictorque THD 
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Figure 14. Line voltage THD 
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7. CONCLUSION 

A SVM-FOC for a three-phase induction motor drive is simulated in MATLAB/Simulink, 
and the dynamic performance is studied too. The starting current, speed regulation, flux distortion and torque 
pulsation are much lower with SVM technique. The simulated responses show that the system performance 
during instant load is good, the simulated torque, current and voltage THD during the steady state 
are improved, with faster dynamic characteristic. In conclusion, SVM-FOC technique has proved good 
performance dynamic for induction motor vector control type. 
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